Magnetophoresis of Flexible DNA-based Dumbbell Structures 
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Controlled movement and manipulation of magnetic micro and nanostructures using magnetic 
forces can give rise to important applications in biomedecine, diagnostics and immunology. We 
report controlled magnetophoresis and stretching, in aqueous solution, of a DNA-based dumbbell 
structure containing magnetic and diamagnetic microspheres. The velocity and stretching of the 
dumbbell were experimentally measured and correlated with a theoretical model based on the forces 
acting on individual magnetic beads or the entire dumbbell structures. The results show that precise 
and predictable manipulation of dumbbell structures is achievable and can potentially be applied 
to immunomagnetic cell separators. 

PACS numbers: 47.63.b,85.70.w,87.14.Gg,47.85.Np 
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Magnetic micro- and nanoparticles find widespread ap- 
plications in biotechnology, e.g. separation of cells and 
biomolecules. Commercially available, micrometer-sized 
magnetic beads (MBs) are polymer spheres containing 
homogenously dispersed superparamagnetic iron oxide 
nanoparticles '1^. To make them suitable for biological 
uses their surface is typically functionalized with high 
affinity binding proteins e.g. Streptavidin or protein A. 
Such affinity functionalized beads can then be bound 
specifically to cells or biomolecules for magnetic sepa- 
ration from biological samples (magnetic pulldown). In 
contrast to such bulk magnetic separation, controlled ma- 
nipulation and movement (magnetophoresis) of individ- 
ual magnetic beads requires more intricate control over 
external magnetic fields. There have been efforts in the 
last 20 years to optimize magnetophoretic conditions for 
this purpose 0, 3, however, applications involving mag- 
netophoresis are still scarce and almost exclusively tested 
on individual MBs [U, i, Q . 

General theoretical basis for MB magnetophoresis can 
be derived by accounting for the forces acting on the MB. 
In the presence of an inhomogeneous magnetic field, a 
magnetic force = V(m • B) is exerted on a MB. Here, 
m is the magnetization of a bead in the magnetic field 
B. In the conditions used in this study the magneti- 
zation of a MB is unsaturated (the saturation magnetic 
field is Bs « 0.5 T), and hence the magnetic response of 
the bead is described as a linear function of the volumet- 
ric magnetic susceptibility, Xb- Therefore, the magnetic 
force becomes: 



V{xb - Xs) 



B(V-B), 



(1) 



where V is MB's volume, Xs is the volume magnetic sus- 
ceptibility of the surrounding aqueous solution and fiQ is 
the vacuum magnetic permeability. In addition to this 
magnetic force, MBs moving in a fluid also experience a 



counteracting viscosity-related drag force. For a spheri- 
cal particle in laminar flow, the drag force is F^ = drnrjv, 
where rj is viscosity of a fluid, and v and r stand for ve- 
locity and radius of a bead, respectively. Thus, the net 
force acting on a MB is given as: 



V{xb - Xs) 



Mo 



(BV ■ B) - GrTTTjv. (2) 
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Magnetic separations of specific immune cell types and 
populations have found widespread use in research, diag- 
nostic and therapeutic uses. Immunochemical methods 
have been developed for attachment of magnetic beads 
to the cell surfaces, typically followed by magnetic sepa- 
ration [6|. Frequent problem encountered in this type of 
"positive" immunomagnetic sorting is that the immune 
cell can "sense" the magnetic bead, thus triggering im- 
munogenic responses such as activation and internaliza- 
tion Q. This often necessitates the use of an alternative, 
"negative" sorting technique that targets for removal all 
cell types with the exception of the cell population of in- 
terest, which remains unattached [1]. Negative sorting is 
an expensive and cumbersome proposition, as it requires 
the use of large quantities and diverse panels of magnetic 
beads directed at many divergent cell types. 

We have reasoned that formation of dumbbell struc- 
tures by attaching magnetic beads to immune cells via 
polymeric linkers would provide sufficient separation be- 
tween the MB and the cell surface to minimize activa- 
tion and internalization, and thus would be an attrac- 
tive alternative to negative immunomagnetic sorting. In 
this letter, we report on an experimental study of con- 
trolled movement of model structures consisting of mag- 
netic beads and latex diamagnetic spheres, joined via 
DNA linkers, as a physical model for magnethophoretic 
separation of immunological cell types via dumbbells. 

An illustration of a dumbbell is shown on Fig. [Ija). 
The magnetic beads M-280 with diameters of 2.8 /im and 
MyOne with diameters of 1 /im (Invitrogen, Oslo, Nor- 
way) were investigated. The iron content for M-280 and 
MyOne was « 12% and w 24%, respectively In this 
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FIG. 1: (a) Schematic illustration of a dumbbell formed from 
a MB, a DNA and a NMB. (b) A disk magnet with a microflu- 
idic channel. The magnetic field is the weakest in the center of 
magnet | B |~ 0.14 T (position x = 5 mm) and the strongest 
I B |~ 0.19 T, at position x = 25 mm. (c) The arm magnet is 
illustrated together with the disk magnet and a microfluidic 
channel. The observation point is roughly indicated. 



study, we concentrated on analyzing the results obtained 
with M-280 MBs. The non-magnetic beads (NMBs) are 
diamagnetic, polymer spheres with a diameter of 5.6 /im 
(Bangs Laboratories, Inc., USA). 

The flexible linker consisted of bacteriophage A-DNA 
(New England Biolabs, USA) with a contour length of 
16.5 /im. A biotinylated oligonucleotide was ligated to 
the first recessed 3' end of A and the purified biotiny- 
lated DNA was then bound to streptavidin coated M-280 
MBs via a strong biotin-streptavidin bond. In order to 
attach a streptavidin-coated NMB to the other terminus 
of the DNA polymer, we used biotin incorporation via 
fill-in of the second recessed 3' end with Klenow poly- 
merase and biotin-dCTP. Finally, the complete DB was 
formed by resuspending the reaction product in aqueous 
buffer and adding streptavidin-coated NMBs. A full de- 
scription of the dumbbell synthesis process will be pub- 
lished elsewhere. Prior to magnetophoresis, the solutions 
were diluted in deionized water containing 1 % by weight 
(1 wt %) of sodium dodecyl sulfate surfactant. This min- 
imizes nonspecific adhesive binding between the beads 
and the surfaces of the channel. 

All experiments on controlled movement and manipu- 
lation of the dumbbells were done in microfiuidic chan- 
nels with widths from 200 /im down to 50 /im. A mould 
with the channel pattern was structured with SU-8 pho- 
toresist on a silicon wafer using standard photolithogra- 



phy. The microfluidic device was cast in polydimethyl- 
siloxane (PDMS) and subsequently sealed on a glass 
cover slip [Toj . 

Magnetophoresis was performed initially on individual 
MBs and subsequently the procedure was applied to the 
dumbbells. The microfluidic chip was loaded with the 
MB suspension (typically 1 /il) and subsequently trans- 
ferred onto a stationary, permanent disk magnet (rare 
earth magnet NdFeB, 50mm diameter, 6mm width). An 
illustration of the magnet with a microfluidic channel is 
shown in Fig. [IJb). The radial magnetic field gradient 
generated by the disk shape (V • B f« 6.5 T/m), enables 
us to achieve directional movement of the beads along 
the channel in a controlled manner (the 2:-direction was 
arbitrarily choosen). We safely neglected effects of the 
gravitational force and thermal energy The move- 
ment of individual MBs was observed with a light mi- 
croscope (51X, Olympus, Japan) equipped with a 20x 
objective (UPlan FT, Olympus, Japan). The trajectories 
were recorded at different spatial positions on the mag- 
net with a digital camera (DP-70, Olympus, Japan). The 
integration time varied between one second (strong mag- 
netic field) up to ten seconds (weak magnetic field) . From 
these images one can calculate the velocity of MBs, by 
dividing the recorded distance with the integration time. 

To estimate the dynamics of a DB in a microfluidic 
channel we considered that the net force acting on a DB 
is composed of the total force given by Eq. [21 acting on 
the MB terminus, while the NMB experiences only the 
drag force. Stretching of the DNA linker occurs only if 
sufficient force difference is exerted on the DB between 
the MB and the NMB. The dumbbell will move with equi- 
librium velocity approximating a value that can be found 
from Eq.Ofor the MB, and which will determine the drag 
force acting on the NMB terminus. The hydrodynamic 
drag force would act as a dynamic anchor, and at equilib- 
rium would be equal to the stretching force. For the disk 
magnet, the stretching force was estimated ~ 0.1 pN pT| . 
which was not sufficient to completely stretch a DB. To 
increase the force, we redesigned our magnetophoretic set 
up, introducing an additional permanent magnet (4.5 mm 
diameter, 2.5 mm width), which can be moved at desired 
position over the stationary disk magnet. A side view of 
a channel sandwiched between the arm and disk magnet 
is shown in Fig.[IJc). This provided a significant increase 
in the magnetic field gradient experienced by a MB. For 
example, even at a distance of 2.5 mm from its center 
(the observation point in Fig. [^c)), the stretching force 
is « 0.5 pN, which should be sufficient to fully stretch a 
single strand of A-DNA 

Figure [U^a) shows an optical image of a microflu- 
idic channel at position x = 10 mm on the disk magnet 
and relatively far away from the arm magnet (> 5 mm). 
Three partially stretched DBs are visible in the channel 
(indicated by circles). The average stretching was found 
to be 1.3±0.3/im. Subsequently, we approached with 
the arm magnet at a distance of « 2.5 mm from the obser- 
vation point. The image shown in Fig.[51[b) is taken at the 
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same position of the microfluidic channel as in Fig. Ufa). 
The movements of all three DBs are clearly visible de- 
spite the short integration time sa 20 msec. Moreover, 
the DBs from left to right along the channel (the arm 
magnet is situated on the right side), show an increase in 
velocity which is reflected in the image as blurring (the 
faster the DB moves, the more pronounced the blurring 
is). The average stretching is found to be 3.1 ± 0.4/im, 
more then doubled compared to the stretching obtained 
from Fig.[2)Ja). This conformal elongation of DB is signif- 
icantly shorter then the full length of A-DNA. The main 
reason for such incomplete streching is the presence of 
multiple DNA molecules in the linker since the large sur- 
face area of the beads contained multiple streptavidin 
molecules. In addition, some of the DNA molecules are 
very likely to be wrapped around the beads, which fur- 
ther reduces the probability of a full contour stretching 
of a DB. It is possible to reduce multiple attachments of 
DNA by sonicating the samples for few seconds (« 4 sec) 
which can lead to shearing of DNA molecules. After 
sonication, the percentage of DBs connected by a sin- 
gle strand was found to be around 20%. However, the 
overall number of DBs was also reduced, and at longer 
sonication times no DBs were visible, due to complete 
shearing of the dumbbell linker. The fluorescent image 
of a DB containing the beads and DNA is shown in inset 
of Fig. Hfb). The image is obtained with a florescence mi- 
croscope where the DNA is stained with a cyanine dimer 
dye, YOYO-1 (Molecular Probes Inc., USA). The numer- 
ous DNA molecules anchored at the edges of the beads 
are clearly visible as a bright halo. The multiple DNA 
binding could be more controllably reduced in the future 
by using magnetic nanoparticles which offer smaller sur- 
face area per particle and hence fewer attachment sites. 

In order to model the system, we needed to know the 
magnetic field and its gradient. The magnetic field was 
measured with a gaussmeter (Lakeshore, 410). From the 
obtained values we were numerically able to extract the 
magnetic field gradient. In Fig. [Hfa), we plotted the ex- 
perimentally obtained velocities at various spatial posi- 
tions on the disk magnet for an individual MB and a 
DB. The velocity-position dependance is simulated by 
numerically solving Eq. [5] using 2.8 /xm as the MB di- 
ameter (2.8 /Ltm and 5.6 /xm for a DB), the viscosity of 
water rj — 8.9 • 10~^ kgm~-'^s~"'^, the magnetic permeabil- 
ity of vacuum fj,o = Att ■ 10~''N/A^, while the magnetic 
susceptibility of the fluid comparative to the MB was ne- 
glected. We would like to emphasize that the volume 
magnetic susceptibility of the MB can be extracted as 
a fitting parameter from the velocity-position measure- 
ments. For an individual MB, we obtain Xb = 0.13 ± 0.04 
which is in good agreement with previously reported val- 
ues [l^ In the case of the DB, the dynamics of 
DNA were neglected, which was justified due to the size 
of the used beads: the drag force associated with the 
beads is greater than for a single DNA molecule, multiple 
DNA molecules in the linker between, or multiple DNAs 
wrapped around beads. From the plot in Fig. [3]^ a), we 
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FIG. 2: (a) An optical microscopy image (20 x objective) 
of a microfluidic channel loaded with DBs at relatively low 
strength magnetic fleld. The three distinct DBs are indicated 
with circles (black, white and dashed lines) . The MBs are po- 
sitioned on the right side of DB due to increase in the gradient 
of magnetic field in the same direction, (b) The same posi- 
tion as in (a) after locaUy increasing the magnetic field and 
the gradient by the arm magnet. Stretched DBs are indicated 
with the same circles as in (a). Inset: A fluorescence image 
(oil immersion 100 x objective, integration time 100 msec), an 
individual, partially stretched DB formed with a DNA as a 
link between the beads. 



see that the experimental results for the MBs and DBs 
are well described by this simple theoretical model. As a 
qualitative check, we also display in Fig.[3]Ja) the product 
I B I (V • B) versus spatial position on the disk magnet 
(dotted line). Clearly, the | B | (V • B) product is the 
main qualitative determinant for the dynamic behavior 
of the objects in this study. Finally, with Fig. Efb) and 
(c) we demonstrate the significant increase in the veloc- 
ity and stretching of a DB, when the arm magnet is used 
in combination with the disk magnet. Fig. [3{b) shows 
histograms of experimentally found velocities, together 
with the corresponding simulation for a single MB, while 
Fig. [3fc) displays similar histograms for the DB. Sim- 
ulated velocities were extracted by calculating the mag- 
netic field and its gradient along the length of the channel 
for the magnetic configuration illustrated in Fig.[ljc) p^ . 
Both histograms in Fig. [3fb) and (c), display a marked 
increase in the velocities of the test objects. While the 
overall trend of velocity increase in the case of the arm 
magnet setup is present, it is in lesser agreement with 
the simulation than in the case of the bare disk mag- 
net. We find that simulation values are highly sensitive 
to the magnetic field gradient and the extracted values 
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strongly depend on the exact choice of observation po- 
sition. Thus, small inherent inaccuracies in precise de- 
termination of the observation point position can lead to 
sizeable disagreement between the model and the data. 

In conclusion, we have investigated magnetophoresis of 
MBs and DBs formed from magnetic and non-magnetic 
beads connected by DNA strands. We find a strong de- 
pendence of dumbbell stretching on the product of the 
magnetic field and its gradient. Experimentally mea- 
sured velocities of the MBs and DBs were compared 
with a theoretical model, assuming simple equilibrium 
between magnetic and drag force. We find solid agree- 
ment between the model and the experimental data. This 
demonstrated approach to use a magnetic force to stretch 
a flexible DNA-based dumbbell - with a hydrodynamic 
drag force acting as a dynamic anchor on the nonmag- 
netic terminus - has the potential to be a versatile tool 
in future applied biological and biomedical separation de- 
vices. 
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FIG. 3: (a) Plot of velocities for the MBs and DBs versus 
their spatial position on the disk magnet. Symbols represent 
measured values, while full lines are simulated velocities cal- 
culated by solving Eq. [2] numerically. The dashed line shows 
a dependance of product | B | (V • B) along the disk magnet 
(right axes), (b) and (c) Histograms of measured and simu- 
lated velocities versus a spatial magnetic position for the arm 
magnet and disk magnet, (b) corresponds to MBs M-280 and 
(c) to DBs. 
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